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ABSTRACT 


The  Navy’s  Impact  Burial  Prediction  Model  creates  a  two-dimensional  time 
history  of  a  bottom  mine  as  it  falls  through  air,  water,  and  sediment.  The  output  of  the 
model  is  the  predicted  burial  depth  of  the  mine  in  the  sediment  in  meters,  as  well  as 
height,  area  and  volume  protruding.  Model  input  consists  of  environmental  parameters 
and  mine  characteristics,  as  well  as  parameters  describing  the  mine’s  release.  The  model 
user  seldom  knows  many  of  these  parameters,  and  those  that  are  known  may  be  of 
questionable  precision. 

In  order  to  determine  which  parameters  had  the  greatest  effect  on  the  model  and 
which  could  be  simplified  or  eliminated,  a  series  of  sensitivity  tests  were  performed.  It 
was  found  that  the  model  data  ingest  could  be  greatly  simplified  without  sacrificing 
accuracy  too  much.  However,  several  parameters  including  sediment  shear  strength  were 
found  to  have  a  large  effect  on  the  model  and  were  investigated  further. 
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L  INTRODUCTION 


In  minesweeping,  success  often  hinges  on  knowing  as  much  as  possible  about  the 
mines  that  have  been  placed  and  the  effects  the  environment  has  had  on  that  placement. 
Since  bottom  mines  cannot  be  searched  for  visually,  and  are  often  difficult  to  locate  with 
conventional  sonar,  an  estimate  of  the  area  or  height  of  the  mine  protruding  from  the 
sediment,  or  the  burial  depth  if  completely  covered,  is  crucial  information  for  the 
planning  and  execution  of  mine  clearance  operations.  Determining  the  likely  mine  burial 
depth  requires  numerical  models  of  the  burial  process  and  knowledge  of  the  environment, 
including  sediment  properties,  waves,  tides,  and  water  depth. 

Ocean  deployed  mines  currently  used  by  the  U.S.  and  other  nations  fall  into  three 
general  categories:  bottom  mines,  moored  mines  and  drifting  mines.  Bottom  mines  rest 
on  the  ocean  floor  and  are  generally  deployed  in  littoral  regions.  Common  placements 
for  bottom  mines  include  shipping  channels,  harbors,  anchorages,  rivers  and  estuaries. 
Bottom  mines  are  deployed  in  one  of  three  ways:  aircraft,  surface  ship  or  submarine. 
Although  mines  are  designed  to  be  deployed  by  a  specific  platform,  most  mines  can  be 
deployed  by  surface  ship  with  little  modification  (NMWEA,  1991). 

Several  numerical  models  have  been  developed  to  simulate  the  mine  burial 
process,  and  constitute  the  only  viable  means  for  determining  a  predicted  burial  depth, 
which  is  critical  information  when  clearing  an  area  of  mines.  The  Impact  Burial  (IB) 
model  was  developed  to  determine  the  depth  at  which  the  mine  comes  to  rest  in  the 
sediment  upon  impact.  Originally  created  by  Amone  and  Bowen  (1980)  at  the  Naval 
Coastal  Systems  Center  in  Panama  City,  Florida,  the  IB  model  was  designed  to  create  a 
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two-dimensional  time  history  of  a  cylindrical  mine  as  it  falls  through  air,  water,  and 
sediment  phases  (Fig.  1).  The  burial  depth  of  the  mine  in  the  marine  sediment  is  then 
calculated  from  the  mine’s  velocity  on  contact  with  the  sediment  and  the  sediment 
characteristics. 

Several  revisions  and  changes  have  been  made  to  the  model  to  refine  the  physics 
and  allow  for  more  realistic  geometry  and  more  extensive  input  from  the  user.  Most 
notable  are  the  changes  made  by  Satkowiak  (1987)  and  Hurst  (1991).  Other  revisions 
involved  translating  to  newer  computer  languages  and  simplifying  the  data  entry  process. 


Figure  1 .  The  trajectory  of  a  cylindrical  mine  as  it  falls  through  three  phases:  air,  water  and 
sediment.  Labels  on  the  right  are  parameters  used  by  the  model  to  calculate  velocity,  attitude, 
and  burial  depth  of  the  cylinder. 
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Currently,  the  model  allows  the  user  to  input  nearly  any  value  for  each 
environmental  parameter.  Many  of  these  parameters  are  rarely  if  ever  known  by  the 
technician,  and  their  inclusion  makes  it  very  difficult  for  the  field  user  to  get  an  accurate 
solution.  With  this  in  mind,  a  sensitivity  test  was  designed  and  executed  with  the 
objective  of  simplifying  the  input  parameters  without  compromising  the  accuracy  of  the 
model’s  output.  By  determining  which  can  be  eliminated  or  simplified  without 
sacrificing  accuracy,  a  complex  model  may  be  made  more  manageable  for  the  customer. 
Several  parameters  were  found  to  lend  themselves  well  to  simplification,  while  others 
had  so  little  effect  on  the  outcome  that  the  default  value  is  sufficient.  Some  parameters, 
however,  have  such  a  large  influence  on  the  final  burial  depth  and  are  known  with  such 
poor  precision  that  a  method  for  more  accurate  field  measurements  should  be 
investigated. 
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II.  DYNAMICAL  ANALYSIS  OF  A 
CYLINDRICAL  MINE 


A.  GEOMETRY  OF  A  CYLINDRICAL  MINE 

As  a  cylinder  falls  through  a  fluid  medium  it  has  a  varying  angle  with  respect  to 

the  vertical  called  its  attitude,  O  (Fig.  2).  The  attitude  of  the  cylinder  determines  the  area 
of  the  cylinder  that  is  perpendicular  to  the  horizontal,  called  the  projected  area.  Ah.  As 
the  attitude  varies  with  time,  the  projected  area  also  changes,  as  do  the  magnitudes  of  the 
vertical  and  horizontal  forces  acting  on  the  cylinder. 


Figure  2.  Horizontal  length,  Lh,  as  a  function  of  cylinder  diameter  and  length,  D 
and  L,  and  attitude  O.  Rotation  rate,  co,  determines  the  attitude  at  a  given  time. 
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B.  ATTITUDE  OF  A  CYLINDRICAL  MINE 

All  forces  and  cavity  calculations  are  affected  by  the  attitude  of  the  cylinder,  O. 

Forces  are  broken  down  into  those  that  act  parallel  to  the  axis  of  the  cylinder,  axial 
forces,  and  those  that  are  perpendicular  to  the  axis,  called  cross  forces.  If  required,  these 
force  values  can  then  be  geometrically  translated  into  horizontal  and  vertical  forces. 

If  the  mass  is  uniformly  distributed  within  a  mine  the  center  of  gravity,  eg,  and 
the  center  of  buoyancy,  cb,  are  the  same.  If  they  differ,  the  moments  about  the  center  of 
buoyancy  which  result  from  forces  in  the  horizontal  and  vertical  directions  are  denoted  as 
Xcg-cb  and  Zcg-cb-  These  values  are  used  along  with  the  weight,  W,  and  velocity,  V,  of  the 
mine  and  the  kinematic  viscosity  and  density  of  the  fluid,  Av  and  p  to  compute  steady 
state  attitude  d>*  for  a  specific  time  increment  (Amone  and  Bowen,  1980) 

0) 

C.  ADDED  MASS 

Following  a  mine's  deployment  by  aircraft,  it  penetrates  the  air,  water,  and 
sediment.  The  periods  during  which  the  mine  falls  through  each  of  these  media  are  called 
the  air  phase,  water  phase,  and  sediment  phase.  A  solid  body  moving  in  a  fluid  originally 
at  rest  behaves  like  a  body  of  increased  inertia  (Mises,  1959),  depicted  using  the  term 
added  mass,  Madded  The  summation  of  the  added  mass  and  the  mass  in  air,  Ma,  is  called 
the  resultant  mass,  Mr: 

Mr  —  Ma  +  Madded  (2) 

Where  Ma  is  the  weight  in  air  (Ma=W/g),  W  is  the  weight  of  the  cylinder  and  g  is  the 
gravitational  acceleration. 
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The  added  mass,  also  called  apparent  mass,  is  a  function  of  the  shape  of  the 
object.  It  is  computed  using  inertial  coefficients  that  account  for  the  flow  processes  about 
the  object  and  the  cylinder  volume,  B.  For  a  cylindrical  object  two  inertial  coefficients 
are  required,  Ki  and  K2,  to  account  for  axial  flow  and  cross  flow: 

Ma  =  K]  p  B  Cos  <E>  +  K2  p  B  Sin  d>  (3) 

D.  PROJECTED  AREA 

The  projected  area  of  a  cylindrical  mine  varies  as  a  function  of  the  cylinder's 
attitude  (Fig.  3).  The  horizontal  and  vertical  lengths  of  a  cylindrical  mine,  Lh  and  Ly,  are 
computed  by: 

Lh  =  L  SinO  +  D  Cos<3> 

Lv  =  L  Cos<E>  +  D  Sind>  (4) 

The  projected  area  for  a  horizontal  cylindrical  mine.  Ah  where  <J>  =  90°,  or  for  a 

vertical  mine,  Av  where  O  =  0°,  is  defined  as: 

Ah  =  LD 

Av  =  B  D2/4 

Thus,  the  projected  area  for  a  cylindrical  mine  with  an  attitude  of  d>  is  computed  by: 

A  =  Ah  Sin2  d>+  Av  Cos2  d>  (5) 
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Figure  3.  Parameters  used  to  compute  projected  lengths  and  projected  area. 

E.  DRAG  COEFFICIENT 

Drag  on  a  circular  cylinder  is  proportional  to  the  sum  of  forces  acting  both 
vertically  and  horizontally.  The  drag  coefficient  used  to  calculate  the  total  drag  must 
therefore  take  each  force  into  account.  The  drag  coefficient,  Cd,  is  derived  from  the 
Reynolds  number  and  the  shape  and  dimensions  of  the  object.  For  a  cylinder,  the  method 
used  depends  on  whether  it  is  a  right  circular  cylinder  or  if  it  has  a  taper  or  rounded  nose. 
Both  cases  are  described  below. 

The  drag  force  on  a  cylinder  is  computed  using  the  projected  area  of  the  immersed 
body.  Ah,  the  velocity  of  the  object,  V,  and  the  drag  coefficient  (Schlichting,  1979): 

Fd  =  0.5  CD  A  V2  A  (6) 
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1.  Drag  Coefficient  for  a  Circular  Cylindrical  Mine 

In  order  to  determine  Cd  we  must  first  determine  the  Reynolds  number  for  the 
object.  A  smooth  right  circular  cylinder  has  a  Reynolds  number  of  less  than  105, 
indicating  that  laminar  flow  is  dominant.  The  Reynolds  number  is  calculated  separately 
for  axial  and  cross  flow  using  the  corresponding  velocity,  V,  and  the  kinematic  viscosity 
of  the  fluid,  Av: 

Re  =  VD/Av  (7) 

The  drag  coefficient  for  the  skin  friction,  Cp,  is  computed  using  the  Reynolds 
number  (Mises,  1959),  which  is  then  used  to  calculate  the  base  drag  coefficient,  C*d: 


Cf  - 

1.328  (Re)-0'5 

(8) 

C*D  = 

0.33D/L  +  Cf  [(3L/D)  +  3(D/L)05] 

(9) 

Further  modification  accounts  for  the  projected  surface  area,  which  will  always  be 
less  than  or  equal  to  the  projected  area  of  a  horizontal  cylinder.  This  is  done  by  applying 
a  ratio  of  Ah  and  the  standard  wetted  area  of  a  cylinder  (equal  to  projected  area  of  a 
horizontal  cylinder),  Aw.  A  generic  coefficient  of  1.1  is  also  applied  to  account  for 
surface  imperfections  on  the  cylinder  (Amone  and  Bowen,  1 980): 

CD  =  CD*  1.1  (Ah/Aw)  (10) 

2.  Drag  Coefficient  for  a  Tapered  or  Round-Nosed  Cylindrical  Mine 

If  the  circular  cylinder  has  a  rounded  nose  or  is  tapered  slightly,  turbulent  flow 

will  be  present  and  the  Reynolds  number  is  computed  by  the  Prandel-Schlichting 

(Schlichting,  1979)  skin-friction  coefficient: 

CF  =  0.455/(log  Re)2  58  (11) 
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The  overall  drag  coefficient  is  calculated  using  both  the  diameter  of  the  base,  Db,  and  the 

maximum  diameter,  D  (Amone  and  Bowen,  1980): 

CD’  =  CF  (Aw/Ah)  *  [1  +  (60(L/D)‘3)  +  0.0025(L/D)] 

CD  =  CD’  + (0.029  (Db/D)3(Cd7°5)  (12) 
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III.  MOMENTUM  BALANCE  OF  A  CYLINDRICAL 
MINE  IN  AIR  AND  WATER 

A.  EQUATION  OF  MOTION 

The  movement  of  a  mine  in  the  air  and  the  water  is  depicted  by  the  momentum 
equation: 

Mr  dV/dt  =  Fw,a  +  Fb  +  Fd  (13) 

Where  V  is  the  velocity  of  the  mine,  Fw>a  is  the  force  due  to  the  air  weight  of  the  mine,  Fb 
is  buoyancy  force  and  Fdis  drag  force  from  equation  (6). 

Buoyancy  force  is  the  upward  force  exerted  upon  a  mine  in  the  gravitational  field 
by  virtue  of  the  density  difference  between  the  mine  and  that  of  the  surrounding  fluid.  We 
use  the  Cartesian  coordinate  system  (x,  y,  z)  with  the  z-axis  in  the  vertical  direction,  and 
use  the  unit  vector  k  along  the  z-axis  (pointing  downward).  The  buoyancy  force  is  then 
computed  using  the  density  value  for  air  or  water,  p,  and  the  volume  of  the  mine,  B: 

Fb  =  -pgBk  (14) 

B.  MINE  ENTRY  INTO  WATER 

A  cylindrical  mine  penetrating  into  water  passes  through  two  distinct  regimes. 

The  first  regime  is  the  cavity  regime.  As  the  mine  pushes  into  the  air-water  interface,  it 
creates  a  cavity  that  consists  of  a  combination  of  air  and  water  particles.  The  ratio  of  air 
to  water  in  the  cavity  decreases  until  the  fluid  properties  become  that  of  water  only,  at 
which  time  the  mine  is  in  the  fully  wetted  regime.  A  temporal  variation  of  the  mine’s 
vertical  position  can  be  calculated  using  a  method  similar  to  that  depicted  in  the  previous 
section. 
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1.  Cavity  Regime 

Upon  impacting  the  water,  the  cylindrical  mine  enters  the  cavity  regime  of  the 
water  phase.  The  fluid  in  the  cavity  is  actually  a  combination  of  air  and  water,  and  the 
properties  change  continuously  with  the  mine's  downward  vertical  motion. 

The  most  critical  properties  of  a  fluid  when  predicting  the  movement  of  an  object 
within  the  fluid  are  density  and  kinematic  viscosity.  In  the  cavity  regime,  the  fluid 
density  is  a  gradient  from  that  of  air,  1.29  kg/m3,  to  that  of  water,  1025  kg/m3.  In  order  to 
describe  the  ratio  of  water  volume  to  total  cavity  volume  we  use  the  term  void  ratio,  fb. 
This  value  is  used  to  determine  the  actual  density  of  the  fluid  in  the  cavity.  Although  it  is 
difficult  to  determine  fD  precisely,  the  gradient  is  assumed  to  be  linear.  When  the  mine  is 
fully  in  the  air  the  void  ratio  is  zero,  and  when  fully  in  water  it  is  one.  The  void  ratio  is 
modeled  based  on  the  assumption  that  the  ratio  of  water  volume  to  air  volume  in  the 
cavity  increases  with  an  increase  in  the  ratio  of  hydrostatic  pressure  to  dynamic  pressure. 
When  the  pressure  generated  by  the  dynamic  loading  of  the  cylindrical  mine  as  it  impacts 
the  water,  called  dynamic  pressure,  is  balanced  by  the  hydrostatic  pressure  of  the  water, 
Ps,  the  cavity  collapses  and  the  mine  is  considered  fully  wetted  (Amone  and  Bowen, 
1980).  The  dynamic  pressure,  Pd,  is  the  summation  of  the  atmospheric  pressure  at  the 
ocean  surface  (assumed  constant)  and  the  pressure  caused  by  the  mine  movement: 

Pd  =  Pa+ViAjFMV2  (15) 

Where  Atfm  is  the  two-phase  mixture  density  in  the  cavity. 

Amone  and  Bowen  (1980)  further  demonstrate  this  relationship  by  stating  that  the 
void  ratio  equals  the  square  root  of  the  pressure  ratio: 

fa  =  JKIK  (16) 
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The  two-phase  mixture  density,  Atfm,  and  kinematic  viscosity,  Atfm,  in  the  cavity 
are  computed  at  each  time  step.  The  fluid  properties  of  density,  pa  and  pw>  and  kinematic 
viscosity,  Aa  and  Aw>  for  air  and  water  are  also  required: 

Atfm  =  Aa  (1  -  fo)  +  foAw  (17) 

Ajfm  —  [(Aa  Aa)  (Aw  Aw)]  /[(l  —  fo)  Aa  Aa  +  fo  Aw  Aw]  (18) 

At  the  each  time  step,  the  hydrostatic  and  dynamic  pressures  are  calculated  from 
the  velocity  of  the  cylindrical  mine,  and  the  void  ratio  is  determined.  When  the 
hydrostatic  and  dynamic  pressures  are  equal,  the  void  ratio  takes  the  value  of  1  and  the 
cavity  collapses.  The  mine  is  then  considered  to  be  in  the  fully  wetted  regime.  At  this 
point  the  velocity  of  the  mine  is  called  the  exit  cavity  velocity.  If  the  actual  water  depth 
is  less  than  the  depth  at  which  the  cavity  collapses,  the  mine  impacts  the  sediment  prior  to 
exiting  the  cavity  and  there  is  no  fully  wetted  regime. 

2.  Fully  Wetted  Regime 

The  trajectory  of  a  mine  in  the  fully  wetted  regime  is  similar  to  that  in  the  air 
except  the  properties  of  the  fluid  have  changed.  The  exit  cavity  velocity  is  used  as  an 
initial  condition  for  the  mine  momentum  equation  (13),  used  to  determine  the  mine 
trajectory.  When  the  vertical  distance  of  the  mine  traveling  in  the  water  equals  the  water 
depth,  the  mine  velocity  is  called  the  bottom  impact  velocity,  which  is  the  initial 
condition  for  determining  the  mine  burial  depth  in  the  sediment. 

C.  MINE  BURIAL  IN  SEDIMENT 

1.  Viscosity  and  Shear  Strength 

Penetration  of  the  cylindrical  mine  into  the  bottom  sediment  depends  primarily  on 
the  attitude  and  velocity  of  the  mine  upon  impact,  as  well  as  the  sediment  properties  of 
density  and  shear  strength.  Initial  impact  of  the  cylindrical  mine  into  the  sediment 
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creates  a  cavity  in  which  the  fluid  properties  of  water  and  sediment  are  interacting.  To 
determine  the  mine  burial  depth  in  the  sediment,  one  must  determine  the  void  ratio  and 
the  two-phase  density  and  kinematic  viscosity.  A  procedure  similar  to  that  depicted  in 
the  previous  section  can  be  used  to  determine  these  parameters.  However,  the 
computation  of  the  two-phase  kinematic  viscosity  for  the  water-sediment  cavity  differs  a 
bit  from  that  of  the  air-water  cavity.  The  kinematic  viscosity  of  the  sediment,  As,  is  not  a 
pure  constant,  but  rather  is  equal  to  the  water  viscosity,  Aw,  plus  that  resulting  from  the 
shear  stress  of  the  sediment: 

As  =  Aw  +  Su  /  (As  d  V/dz)  (19) 

Where  As  is  the  density  of  the  sediment  and  Su  is  the  shear  strength. 

2.  Mine  Burial  Dynamics 

The  vertical  momentum  balance  of  a  mine  in  the  sediment  phase  is  given  by: 

Mr  dV/dt  =  Fw,a  +  Fb  +  Fd  +  Fc  +  F,  (20) 

where  Fb  is  the  buoyancy  force  in  the  sediment,  Fc  is  the  compressive  force,  and  Fs  is  the 
shear  force.  Fc  and  F,  are  additional  forces  (different  from  air  and  water  phases)  exerted 
on  the  mine  by  the  sediment.  They  are  proportional  to  shear  strength  of  the  sediment  and 
the  projected  area  of  the  mine.  If  the  mine  is  a  right  circular  cylinder,  the  compressive 
force  is  twice  the  shear  force: 

FC  =  2FS 

Fs  =  SuA  (21) 

The  mine  burial  depth  is  predicted  by  integrating  (20)  with  respect  to  time  until 
the  mine  velocity  becomes  zero.  Accurate  values  for  sediment  properties  are  essential  to 
the  accuracy  of  this  process.  Shear  strength  and  density  have  a  strong  impact  on  the 
computation  of  all  forces  as  well  as  buoyancy  weight  and  added  mass. 
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IV.  MINE  IMPACT  BURIAL  PREDICTION  MODEL 


A.  DEVELOPMENT  HISTORY 

Based  on  the  fluid  dynamics  of  a  cylindrical  mine,  Amone  and  Bowen  (1980) 

developed  a  mine  impact  burial  prediction  model  called  the  Impact  Burial  Prediction  (IB) 
Model.  The  basic  model  created  a  two-dimensional  free-fall  history  of  a  right  circular 
cylinder  as  is  fell  through  air,  water  and  sediment  phases  and  calculated  the  burial  depth 
of  the  mine  after  it  came  to  rest  in  the  sediment.  The  initial  model  had  a  number  of 
shortcomings,  particularly  when  handling  soft  or  hard  sediments  or  certain  environmental 
conditions.  It  also  handled  the  water-sediment  cavity  regime  poorly. 

Satkowiak  (1987)  modified  the  IB  model,  including: 

■  Corrected  reference  flow  used  in  drag  calculations 

■  Corrected  added  mass  term  equation 

■  Reworked  equations  for  sediment-cavity  regime 

■  Allowing  for  drag  due  to  cylinder’s  nose 

■  Allowing  for  rounded  noses 

■  Inclusion  of  water  temperature 

■  Inclusion  of  retarding  forces  in  semi-solid  sediment 


The  IB  model  was  extensively  revised  by  Hurst  (1991).  Equations  for  forces 
acting  on  the  mine  were  redefined  and  equations  to  simulate  rotational  movement  of  the 
mine  were  added.  Five  primary  areas  were  addressed  and  refined: 

■  Calculation  of  fluid  drag 

■  Calculation  of  forces  at  air-water  interface 

■  Calculation  of  forces  upon  impact  with  the  sediment 

■  Allowing  for  rotational  movement  of  mine  as  it  falls 

■  Allowing  for  multilayered  sediments  (formerly  only  the  deepest  point  was 
considered) 

Also  added  in  1991  was  the  ability  of  the  model  to  calculate  fall  angle  dynamically,  as  it 
would  in  reality. 


15 


In  1993,  Mulheam’s  formulation  for  sediment  bearing  strength  (1993),  which 
takes  sediment  shear  strength  as  well  as  the  object’s  mass  and  geometry  into  account,  was 
added  to  the  IB  model.  The  model  was  also  modified  at  this  time  to  allow  for  displaced 
centers  of  mass. 

Currently,  the  twenty-sixth  edition  of  the  IB  model  is  available  in  ANSI  C  format. 
The  model  features  interactive  data  input  with  the  user  being  prompted  at  every  step.  No 
user  manual  is  required,  and  none  exists  at  this  time. 

B.  IB  MODEL  INPUT  PARAMETERS 

The  EB  model  has  a  variety  of  parameters  that  must  be  entered  by  the  user  prior  to 
running  the  model.  Several  of  these  parameters  are  unknowns,  even  to  the  mine 
sweepers,  and  others  are  easily  determined  based  on  intelligence  and  published  charts  or 
atlases. 

One  possible  unknown  parameter  that  is  required  is  mine  type.  If  the  specific 
type  of  mine  is  known,  the  user  enters  the  mine’s  characteristics  of  air  weight,  water 
weight,  length,  maximum  diameter,  taper  length,  base  diameter,  and  center  of  mass  offset 
distance  from  the  published  mine  lists.  If  unknown,  an  educated  guess  must  be  made  and 
a  profile  selected  from  those  available  in  the  model  (Table  1).  Most  of  the  mine  profiles 
available  in  the  model,  and  most  mines  used  today,  are  right  cylinders  with  a  center  of 
mass  in  the  geometric  center  of  the  mine.  The  model  also  has  the  capability  to  handle 
tapered  mines  and  offset  centers  of  mass.  Both  the  weight  and  geometry  of  the  mine  have 
a  large  effect  on  the  model’s  output,  so  care  should  be  taken  in  selecting  the  mine  type 
(Fig-  4). 
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Mine  Name 

Mine 

Case 

Number 

■a 

m 

'./■'[iiiSjt  -i 

Base 

Diameter 

(m2) 

Mine  A 

1 

612 

354 

1.524 

0.48 

0 

0.48 

MineB 

2 

227 

170 

1.676 

0.292 

0.396 

0.076 

Mine  B-l 

3 

227 

0.292 

0.396 

0.076 

Mine  B  mod 

0 

4 

499 

181 

1.767 

0.47 

0 

0.47 

MineC 

5 

862 

408 

2.134 

0.559 

0 

0.559 

Mine  A  mod 

1 

6 

442 

1.524 

0.48 

0 

0.48 

PSIMine 

7 

515 

1.481 

0.48 

0 

0.48 

Korean  Mine 

8 

538 

251 

1.49 

0 

Bowen  Mine 

9 

964 

413 

2.316 

0.597 

0 

0.597 

NZ  Dummy 
Mine 

10 

457 

263 

1.219 

0.445 

0 

US 

Penetrometer 

11 

3.36 

0  • 

Table  1 :  Parameters  for  mines  that  are  available  for  selection  in  the  model.  Center  of  mass  offset 
difference  can  also  be  specified  (in  meters)  but  is  0  for  all  selections  available  in  the  IB  model. 
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The  initial  orientation  of  the  mine  can  be  set  between  0°  and  90°,  with  0°  being  a 
vertical  initial  orientation  and  90°  being  horizontal.  This  parameter  may  be  based  on  the 
knowledge  of  how  the  mines  were  laid,  but  if  unknown  the  default  of  90°  will  yield  the 
minimum  burial  depth.  Another  unknown  value  is  the  rotation  rate,  the  constant  rate  a 
mine  would  rotate  in  degrees  per  second.  Because  rotation  rate  is  known  only  in 
theoretical  cases,  the  default  has  been  set  at  zero. 

If  dropped  vertically  with  no  initial  rotation  rate,  the  mine  will  still  tend  to  acquire 
a  horizontal  velocity  component  as  it  falls  through  the  air  and  water  phases.  This  results 
in  the  mine  impacting  the  sediment  at  a  random  angle,  significantly  altering  the  actual 
burial  depth.  Satkowiak  (1988)  suggested  that  since  the  attitude  of  the  mine  upon  impact 
with  the  sediment  cannot  possibly  be  known,  even  for  controlled  tests,  the  model  should 
be  used  mainly  to  determine  the  maximum  and  minimum  limits  of  the  burial  depth.  Since 
0°  initial  orientation  provides  the  shallowest  burial  and  0°  the  deepest,  the  limits  can  be 
easily  determined  and  the  range  of  burial  depth  used  with  confidence. 

C.  AIR  AND  WATER  PHASES 

As  the  mine  passes  through  the  air  and  water  phases,  all  forces  acting  upon  it  are 
calculated  and  summed.  The  acceleration  is  then  calculated  and  integrated  to  provide  the 
velocity  of  the  mine  as  it  enters  into  the  next  phase.  The  calculations  for  the  air  and 
water  phases  are  similar,  and  for  the  most  part  are  handled  by  a  generic  fluid  subroutine. 
The  first  thing  computed  at  each  time  step  by  the  model  is  the  effective  axial  and  cross 
masses  in  air,  cross  and  axial  velocity,  and  the  distance  from  the  bottom  of  the  mine  to 
the  center  of  mass  and  center  of  gravity  for  the  current  attitude  of  the  mine.  These  values 
are  then  used  to  calculate  the  drag  coefficients,  drag  force  and  torque  on  the  mine. 
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Finally,  the  velocities  are  translated  to  x  and  z  movement  (in  the  Cartesian 
coordinate  system)  and  O  is  recomputed  to  represent  <$*,  the  steady  state  attitude  at  the 
end  of  the  time  increment.  All  values  are  calculated  for  a  time  increment  of  0.01 
seconds,  and  the  values  are  printed  to  the  screen  in  increments  of  0.5  seconds. 

The  interface  between  each  phase  has  a  cavity  regime  where  the  nose  of  the  mine  is  in  a 
cavity  that  it  actually  a  mixture  of  the  two  phases.  The  composition  of  the  mixture 
changes  rapidly  as  the  mine  pushes  through  the  cavity.  Depth  of  the  cavity  is  calculated 
from  the  previous  time  step’s  values  for  attitude  and  depth  of  mine.  The  cavity  exists 
until  dynamic  pressure  equals  hydrostatic  pressure,  at  which  point  the  cavity  collapses. 

D.  SEDIMENT  PHASE 

The  physical  properties  of  the  sediment  used  by  the  model  are  the  depth,  density 
and  shear  strength  of  each  layer.  This  information  may  be  available  for  some  areas,  but 
no  database  currently  exists  that  contains  sediment  profiles  including  shear  strength.  This 
value  is  independent  of  density  in  saturated  sediment  and  is  difficult  to  measure, 
especially  for  soft  sediments  where  a  sample  profile  is  extremely  difficult  to  obtain. 
There  are  three  sediment  profiles  available  in  the  model  that  represent  hard,  medium  and 
soft  sediment  types.  All  mines  and  all  initial  conditions  produce  complete  burial  in  the 
soft  sediment  case.  The  most  variance  in  burial  depth  comes  from  the  medium-density 
sediment  profile,  and  for  this  reason  it  was  used  in  this  sensitivity  study  unless  otherwise 
noted. 

The  main  contributions  on  the  slowing  of  the  mine  as  it  impacts  the  sediment, 
according  to  Hurst,  are  bearing  strength  of  the  sediment  (70%),  hydrodynamic  drag  of  the 
sediment  (25%),  and  buoyancy  in  the  sediment  (5%). 
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1.  Bearing  Strength 

In  the  original  IB  model,  bearing  strength  was  assumed  to  be  the  undrained  shear 
strength  times  ten.  It  was  soon  determined  that  this  simple  assumption  caused  the  burial 
depth  to  be  consistently  underestimated.  By  using  a  method  devised  by  Mulheam  (1993) 
for  calculating  the  bearing  strength  that  incorporates  the  mine  dimensions  and  the  impact 
velocity,  a  more  accurate  result  is  obtained.  Mulheam  bearing  strength,  BSM,  is 
calculated  using  sediment  shear  strength,  Su,  and  the  rate  of  loading  factor,  rlf: 

rlf  =  (0.5  (mine  impact  velocity)  /  D)° 15 

BSm  =  5. 14  Su  (1  +  D  /  (5. 14Lh))  ( 1  +  0.4z  /  D)  (rlf)  (22) 

Where  Lh  is  the  horizontal  projection  of  mine,  D  is  the  mine’s  diameter  and  z  the  depth 
from  the  mine’s  lower  surface  to  the  water-sediment  interface. 

In  the  IB  model,  this  formulation  can  be  toggled  on  or  off  by  the  user.  Best 
agreement  between  measured  and  modeled  burial  depth  was  found  when  this  formulation 
was  used  and  shear  strength  was  assumed  to  be  constant  for  depths  greater  than  0.55 
meters.  The  correct  formulation  for  the  rate  of  loading  factor  is  still  uncertain,  however. ' 

2.  Hydrodynamic  Drag  in  Sediment 

Hydrodynamic  drag  in  the  sediment  is  the  force  required  to  push  the  sediment 
aside.  Assuming  a  low  Reynold’s  number  of  105,  accurate  for  a  right  circular  cylinder, 
and  using  a  standard  cavitation  correction  of  0.55,  a  drag  coefficient  of  0.98  can  be 
empirically  derived.  Cdh,  the  hydrodynamic  drag  coefficient,  is  then  calculated  for  the 
vertical  velocity,  Cz,  at  each  time  increment: 

Cdh  =  (0.98  *Cz2)/2  (23) 

Cdh  and  the  calculated  area  in  contact  with  the  sediment  are  then  used  to  calculate  the 
force  this  drag  has  on  the  cylinder. 


20 


3.  Buoyancy  in  Sediment 

Buoyancy  in  the  sediment  is  due  to  the  cavity  formed  as  the  nose  pushes  against 
the  sediment  and  forms  a  cavity.  Buoyancy  is  a  result  of  hydrostatic  pressure  and  is 
calculated  using  the  density  of  the  sediment  ps,  density  of  the  mine  pm,  and  the  depth  in 
the  sediment  h: 

Fb  =  g  h  (ps  -  Pm)  (24) 

It  is  assumed  that  the  hydrostatic  forces  work  uniformly  on  the  surface  area  of  the  leading 
edge  of  the  mine  that  is  in  contact  with  the  sediment.  The  total  buoyancy  force  is  then 
calculated  using  the  surface  area  of  the  leading  edge  and  the  buoyancy  coefficient  above. 

E.  ENVIRONMENTAL  PARAMETERS  IN  THE  IB  MODEL 

The  altitude  from  which  the  mine  is  released  determines  the  velocity  and  attitude 

of  the  mine  as  it  reaches  the  air-water  interface.  If  a  mine  does  not  fall  straight  down  but 
rather  “tumbles”  with  a  constant  rate  of  rotation,  simulated  in  the  model  by  providing  a 
rotation  rate  0,  the  attitude  of  the  mine  upon  reaching  the  water  is  impacted  greatly  by  the 
release  altitude.  Although  not  accounted  for  in  the  model,  this  rotation  rate  may  be 
caused  or  affected  by  wind. 

In  the  water  phase,  this  rotation  rate  is  damped  significantly.  However,  it  still  has 
a  great  effect  on  the  angle  the  mine  makes  with  the  sediment  upon  impact.  Currents  may 
affect  the  rotation  rate  in  the  model,  but  again  are  not  accounted  for  in  the  model.  The 
water  depth  only  has  an  effect  on  impact  velocity  if  it  is  less  than  that  required  for  the 
mine  to  reach  terminal  velocity,  the  velocity  at  which  the  deceleration  due  to  frictional 
drag  is  equal  to  the  acceleration  from  gravity.  The  velocity  at  which  this  equilibrium  is 
reached  is  a  function  of  the  weight  of  the  mine.  Since  mines  are  laid  in  shipping  channels 
almost  exclusively,  one  may  assume  that  water  depths  in  excess  of  that  required  for  a 
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mine  to  reach  terminal  velocity  are  the  norm.  Water  temperature  has  an  effect  on  the 
viscosity  of  seawater,  and  hence  increases  the  drag  of  the  seawater  on  the  mine. 

Properties  of  the  sediment  are  represented  by  density  and  shear  strength  profiles. 
Density  of  marine  sediment  tends  to  have  a  s-shaped  profile  with  sharper  gradients  as 
density  increases.  Shear  strength,  the  ability  of  the  sediment  to  withstand  pressure 
without  deforming,  also  typically  has  a  s-shaped  profile  and  increases  with  distance  from 
the  water  interface.  The  shear  strength  is  related  to  the  level  of  cohesion  between  the 
sediment  particles.  The  density  range  of  concern  to  the  mine  impact  burial  problem  is 
1375  to  1600  kg/m  .  Factors  contributing  to  shear  strength  are  the  type  of  material,  water 
content,  history  of  stress  or  disturbance  and  time  since  deposition  (Noorany,  1985). 
Although  both  increase  with  distance  from  the  interface,  there  is  no  clear  correlation 
between  shear  strength  and  density.  Figure  5  is  a  scatter  plot  of  density  and  shear 
strength  values  for  62  sediment  samples,  all  taken  at  the  water-sediment  interface.  For 
this  particular  data  set,  the  correlation  is  extremely  weak.  Shear  strength  at  the  water- 
sediment  interface  can  be  measured  in  situ  with  a  vane  penetrometer  or  other  instrument. 

A  profile  of  shear  strength  such  as  is  called  for  in  the  model  must  be  measured 
from  a  core  sample  in  the  laboratory.  This  process  is  time  consuming  and  expensive,  and 
no  database  of  shear  strength  values  currently  exists.  The  term  bearing  strength,  as  used 
in  the  IB  model  data  ingest,  refers  to  the  undrained  shear  strength  times  10.  This  value, 
however,  is  converted  back  to  shear  strength  and  used  in  the  Mulheam  bearing  strength 
equation. 


22 


Although  not  a  parameter  considered  in  the  IB  model,  wave  action  has  a  direct 
effect  on  water  depth  and,  therefore,  on  velocity  of  the  mine  as  it  reaches  the  sediment 
interface.  This  effect  only  becomes  significant  when  the  ratio  of  water  depth  to  wave 
height  is  high,  and  only  at  very  low  release  altitudes  (Table  2). 


Wave  Height 

Mean  Water  Depth 

1  meter 

2  meter 

5  meter 

5  meters 

0.12 

0.21 

13.7  meters 

0.03 

0.06 

0.16 

20  meters 

0 

0 

0 

Release  altitude=l  .5m 

Wave  Height 

Mean  Water  Depth 

1  meter 

2  meter 

5  meter 

5  meters 

0.02 

0.05 

13.7  meters 

0 

0 

0 

20  meters 

0 

0 

0 

Release  altitude=150m 

Table  2.  Effect  of  wave  height  on  burial  depth.  Wave  heights  were  set  around  a  mean  water 
depth  and  release  altitude  was  selected  to  represent  a  ship  or  aircraft  delivery.  Tf  released  from 
150  meters,  wave  height  has  little  effect.  If  released  from  near  the  water  surface,  wave  height 
upon  release  has  a  significant  impact  if  the  depth/wave  height  ration  is  high. 
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V.  SENSITIVITY  STUDY  FOR  THE 
CYLINDRICAL  MINE 

The  purpose  of  this  study  is  to  ascertain  which  variables  the  model  is  most 
sensitive  to  and  which  can  be  simplified  or  eliminated  in  order  to  simplify  its  use.  Since 
some  variables  are  typically  unknown  by  the  user,  it  is  important  to  determine  which  of 
these  have  the  most  impact  on  the  model  and  which  can  be  reduced  to  toggled  values  or 
default  values  without  greatly  impacting  the  accuracy  of  the  model.  The  model  was 
altered  to  allow  most  parameters  to  be  set  and  a  loop  run  of  one  variable  at  a  time.  The 
range  of  each  variable  was  set  to  represent  all  possible  conditions  the  model  would  be 
used  under.  It  should  be  noted  that  wind  and  currents  are  not  accounted  for  in  this  model. 
However,  the  only  impact  they  would  have  would  be  on  the  attitude  of  the  mine  as  it 
enters  each  phase.  All  runs  were  made  with  preset  mine  profile  “Korean  Mine”,  which 
has  a  dry  weight  of  538  kg,  a  wet  weight  of  251  kg,  and  a  uniform  diameter  of  .475 
meters. 

Since  the  model  calculates  burial  depth  and  then  geometrically  calculates  the  height, 
volume  and  area  protruding,  these  values  are  proportional.  To  confirm  this,  we  created 
derivative  plots  of  these  values  for  one  case  and  found  the  shapes  of  the  curves  to  be  very 
similar  (Fig.  6).  Burial  depth  is  used  to  explain  most  of  the  sensitivity  test,  except  where 
height  protruding  is  more  descriptive. 

A.  SENSITIVITY  TO  RELEASE  MEDIUM  PARAMETERS 

Figure  7  demonstrates  the  variation  of  the  release  medium  parameters  of  altitude, 
water  depth  and  water  temperature.  Altitude,  when  varied  from  0  to  1000  meters,  has  a 
small  impact  on  burial  depth  (relative  difference  of  18%).  When  a  more  realistic  upper 
limit  of  300  meters  for  a  mine  laying  aircraft  is  applied,  the  relative  difference  drops  to 
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just  9%.  Water  depth  has  an  effect  on  the  burial  depth  only  if  less  than  the  depth  needed 
for  a  mine  to  reach  terminal  velocity,  in  this  case  20  meters.  At  depths  greater  than  this 
value,  the  mine  reaches  terminal  velocity  in  the  seawater  and  excess  water  depth  has  no 
effect  on  burial  depth.  At  depths  from  0  to  20  meters,  the  variance  in  burial  depth 
depends  on  both  altitude  and  water  depth  since  the  vertical  velocity  of  the  mine  as  it 
enters  the  water  becomes  pertinent  (Fig.  8). 

Water  temperature  was  found  to  have  no  effect  on  the  model’s  outcome. 
Although  temperature  variance  does  alter  the  density  of  water  up  to  3%  and  also  affects 
the  viscosity  (Stanley,  1969),  this  effect  is  not  significant  enough  to  alter  the  burial  depth 
value  calculated  by  the  model. 

All  cases  discussed  thus  far  assumed  an  initial  angle  of  90°  with  respect  to  the 
horizontal  and  with  a  rotation  rate  of  zero.  This  produces  a  situation  where  the  mine  is 
heading  directly  downward  throughout  the  entire  simulation,  resulting  in  the  maximum 
burial  depth.  When  this  initial  attitude  is  varied,  the  burial  depth  is  affected  greatly  as 
outlined  in  Table  3. 


Altitude  =  1.5  meters 

Altitude  =150  meters 

Fall  Angle  =  0° 

0.977  m 

2.405  m 

Fall  Angle  =  90° 

0.342  m 

0.359m 

Table  3.  Maximum  and  minimum  burial  values  for  a  mine  released  from  1.5  or  150 
meters.  An  initial  fall  angle  and  subsequent  sediment  impact  angle  of  0°  indicates  a 
perpendicular  orientation  and  maximum  burial  depth.  Fall  angle  of  90°  indicates  the 
mine  is  parallel  to  the  sediment  and  yields  a  minimum  burial  depth. 
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B.  SENSITIVITY  TO  SEDIMENT  CHARACTERISTICS 

As  expected,  sediment  parameters  are  the  most  critical  element  in  determining 

how  deep  the  mine  was  buried  when  it  came  to  rest.  Sensitivity  to  the  alteration  of 
sediment  density  and  shear  strength  was  tested  two  ways.  First,  six  sediment  profiles 
were  entered  into  the  model  and  the  resulting  burial  depth  was  examined  (Fig.  9).  These 
included  three  profiles  from  Sydney  Harbor  (Mulheam,  1993)  and  three  profiles  available 
for  selection  in  the  IB  model.  The  profiles  included  in  the  model  are  called  simply 
“softsed”,  “medsed”,  and  “hardsed”  and  do  not  clearly  correspond  to  specific  sediment 
types.  Second,  simplified  cases  of  a  single  layer  of  sediment  were  used  with  constant 
density,  varying  shear  strength  and  constant  shear  strength  with  varying  density. 

1.  Sensitivity  to  Shear  Strength  and  Density 

Figure  10  illustrates  the  sensitivity  of  burial  depth  on  density  and  shear  strength. 
Here,  a  simple  profile  of  just  one  layer  was  used  and  density  and  shear  strength  were 
varied  separately.  All  other  parameters  were  kept  unchanged  as  default  values.  Plot  (a)  is 
the  burial  depth  with  shear  strength  held  constant  and  varying  density  from  1000  to  2000 
kg/m3.  Shear  strength  of  1  kPa  indicates  extremely  soft  sediment,  and  density  has  a 
noticeable  effect  on  burial  depth  of  37%.  At  the  more  common  shear  strength  values  of  5 
to  15  kPa  density  has  little  effect,  just  3.7%.  Plot  (b)  illustrates  the  effect  of  varying 
shear  strength  while  keeping  density  constant.  Again,  we  see  the  greatest  impact  of 
density  value  on  the  model  output  at  low  shear  strength  values.  As  shear  strength 
increases,  so  does  the  influence  of  varying  density. 

Tables  5  and  6  illustrate  the  effect  on  burial  depth  of  simplifying  the  sediment 
profiles  in  two  ways.  Several  methods  for  simplifying  the  sediment  profile  requirements 
were  investigated,  using  the  full  profile  case  as  a  control.  First,  the  density  and  shear 
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strength  were  held  constant  to  5  meters.  The  relative  difference  is  under  26%  for  all 
profiles  using  this  simplification.  Next,  a  process  of  manufacturing  sediment  profiles 
using  the  density  values  measured  at  the  water-sediment  interface  was  derived  and 
applied  to  the  model.  The  profiles  were  assumed  to  consist  of  a  constant  value  layer  at 
the  surface  to  a  depth  hi,  a  sharp  gradient  to  h2,  and  then  constant  to  a  depth  of  5  meters. 
The  profiles  were  first  applied  to  density  only,  holding  shear  strength  constant,  and  then 
to  both  density  and  shear  strength.  Values  for  hi,  h2,  p(h2),  and  x(h2)  were  calculated  by 
applying  p0  and  xa  to  polynomials  derived  from  the  data.  The  softsed  and  medsed 
profiles  create  the  greatest  differences  from  the  control  in  all  cases. 

Interestingly,  creating  a  simulated  density  profile  and  keeping  the  shear  strength 
value  constant  had  no  effect  on  the  burial  depth  result  when  compared  to  keeping  both 
values  constant  for  five  meters.  This  serves  to  underscore  the  fact  that  it  is  the  shear 
strength  parameter  that  has  the  primary  influence  on  burial  depth,  not  the  more  easily 
measured  density  parameter. 

2.  Simulated  Sediment  Profiles 

Several  attempts  were  made  to  manufacture  shear  strength  profiles  from  density 
and  shear  strength  values  measured  at  the  interface.  This  was  explored  in  order  to 
determine  if  a  viable  method  of  simplifying  the  data  entry  for  the  sediment  phase  could 
be  devised.  One  attempt  consisted  of  applying  a  fitted  polynomial  to  measured  density 
and  shear  strength  values  to  create  a  synthetic  profile  from  only  interface  values.  Values 
for  the  sediment  profiles  used  in  the  study,  calculated  with  the  following  equations,  are 
listed  in  table  5. 
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Based  on  the  density  profiles  in  figure  9,  we  empirically  derived  curve-fitting 

equations  to  represent  the  density  profile, 

hi  =  -0.000061833*po+  0.01609 

h2  =  -0.0015*po  +  3.10 

Su(h2)  =  [(2.2048*  10'5)po4]  -  (0.109p3)  +  201.1381p2[(1.6432*105)p]  + 

(5.0136*107) 
p(h2)  =  po  +  250 


1 

Rose  Bay 

Off 

Woolwich 

0.083 

0.05 

0.084 

1.21 

0.40 

1.23 

I  J0558MI 

1260 

1800 

1250 

1510 

2050 

1500 

■m i 

9 

13 

8 

t(h2) 

44 

20.9 

41 

softsed.sed 

medsed.sed 

hardsed.sed 

MNOSSHK 

0.093 

0.084 

0,081 

1.45 

1.23 

1.15 

SS5SSB 

1100 

1250 

1300 

_ p(k) _ 

1350 

1550 

HBHH 

1 

1.5 

13 

t(h2) 

3.5 

10 

20.9 

Table  4.  Values  calculated  using  equations  derived  by  fitting  a  polynomial  to  known 
density  and  shear  strength  profiles.  The  precision  of  this  simplified  method  of  depicting 
a  sediment  profile  is  demonstrated  in  table  5. 


31 


Mossman 

Bay 

Rose  Bay 

Off 

Woolwich 

Full  Profile 

-- . -  . --- . 

0.103  m 

0.074  m 

0.115m 

Measured  pto  S„  held 
constant  to  5  meters 

0.121m 

0.093  m 

0.132  m 

Relative  Difference 

15% 

20% 

13% 

Density  and  Shear 
strength  profiles 
created  using 
measured  p0  and  Suo 

0.101m 

0.059  m 

0.105  m 

Relative  Difference 

2% 

20% 

1% 

softsedsed 

medsed.sed 

hardsed.sed 

Full  Profile 

0.523  m 

0.342  m 

0.084  m 

Measured  Po,  Su  held 
constant  to  5  meters 

0.683  m 

0.463  m 

0.094  m 

Relative  Difference 

23% 

26% 

11% 

Density  and  Shear 
strength  profiles 
created  using 
measured  p0  and  Suo 

0.300  m 

0.179  m 

0.085  m 

Relative  Difference 

43% 

48% 

i% 

Table  5.  Predicted  burial  depths  using  manufactured  profiles  based  on  measured 


values  at  the  interface. 
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Hayter  (1986)  discussed  an  equation  originally  derived  by  Krone  (1963)  for 
deriving  shear  strength,  Su,  from  density  using  empirically  derived  coefficients  a  and  (5: 

Su=app 

Values  for  a  and  (3  must  be  calculated  for  each  separate  sediment  type,  after  which  the 
shear  strength  can  simply  be  calculated  using  the  coefficients.  Figure  1 1  illustrates  the 
impact  of  varying  a  and  P  on  the  model  output,  given  a  constant  density.  The  profile  was 
assumed  to  consist  of  one  layer  of  homogenous  material.  As  expected,  as  a  and  P 
increase,  shear  strength  also  increases  and  burial  depth  decreases.  Figure  12  is  a  series  of 
contour  plots  with  varied  values  as  the  axes.  The  contours  represent  predicted  burial 
depth  values.  Known  shear  strength  values  are  marked  on  the  corresponding  density  plot. 
For  all  cases,  there  is  a  unique  value  of  the  coefficients  that  will  produce  a  shear  strength 
value  given  a  specific  density.  Please  note  that,  while  they  are  plotted  here  as  one  density 
value  and  one  shear  strength  value  per  sediment  type,  a  change  in  density  would  produce 
a  corresponding  change  in  shear  strength  that  could  be  determined  by  use  of  the  same  two 
unique  values  of  the  a  and  P  coefficients. 
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Derivatives  of  burial  depth  and  height  protruding 


(a) 


Figure  6.  Derivative  plot  of  all  four  output  parameters  from  one  IB 
model  run.  Notice  the  shape  of  each  derivative  plot  is  similar, 
confirming  that  one  can  be  used  to  represent  sensitivity  levels  of  all 
output  parameters.  Burial  depth  was  examined  in  the  study  for  this 
paper,  as  it  is  the  primary  result  obtained  by  the  model. 
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Derivative  of  volume  exposed  Derivative  of  area  exposed 
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Water  Depth  (m) 

Figure  6.  Continued  from  previous  page. 
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Figure  7.  Effect  of  (a)  release  altitude  (m),  (b)  water 
depth  (m)  and  (c)  water  temperature  (°C)  on  predicted 
burial  depth  (m).  Values  were  preliminarily  chosen  to 
represent  all  conditions  under  which  the  IB  model 
may  be  used. 
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Burial  Depth  (m) 


Water  Depth  (m)  Release  Altitude  (m) 

Figure  8.  Three-dimensional  plot  of  (a)  burial  depth  (m)  and  (b)  height 
protruding  (m)  as  both  release  altitude  (m)  and  water  depth  (m)  are 
varied.  Height  protruded  is  illustrated  here  to  clarify  the  levels  at  which 
these  parameters  become  less  influential  in  the  IB  prediction. 
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Depth  (m) 
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Figure  9.  Sediment  profiles  of  (a)  density  (kg/m3)  and  (b)  shear  strength 
(kPa)  used  in  the  sensitivity  study.  Data  obtained  from  Mulhearn 
(1993). 
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Shear  Strength  (kPa) 

Figure  9.  Continued  from  previous  page. 
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Burial  Depth  (m) 


Burial  Depth  (m)  Burial  Depth  (m) 


s 


alpha=5.4*l0 


alpha=6. 75*10 


Figure  11.  Effect  of  a  and  (3  coefficients  on  predicted  burial  depth  (m), 


beta  beta 


rho=1000 


rho=1 250 


Figure  12.  Contour  plots  of  predicted  burial  depth  (m)  with 
respect  to  a  and  p  coefficients.  Asterisks  mark  the  values  of  a 
and  P  that  correspond  to  the  interface  values  of  the  three 
sediment  profiles  from  Mulhearn  (1993)  used  in  this  study  and 
to  a  subset  of  the  Voelkner  (1973)  data  set. 
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VI.  RECOMMENDATIONS 


As  expected,  there  are  several  parameters  that  are  both  rarely  known  by  the 
operator  and  of  little  import  to  the  outcome  of  the  model.  In  order  to  make  the  model 
easier  to  Use  with  out  sacrificing  accuracy,  these  parameters  should  be  simplified  as  much 
as  possible. 

Water  temperature  was  found  to  have  no  effect  on  burial  depth,  and  should  be 
eliminated  from  the  list  of  variables.  Altitude  values  should  be  simplified  to  represent  a 
mine  laying  platform  and  the  most  likely  height  for  release  from  that  platform.  For 
instance,  it  is  more  likely  that  an  operator  would  know  if  the  mines  were  laid  by  a  ship  or 
by  a  certain  type  of  aircraft  that  the  enemy  has  than  at  which  altitude  that  aircraft  was 
flying.  Using  intelligence  and  e.xperience,  one  could  form  a  rough  database  for  selected 
platforms.  The  difference  between  an  altitude  of  1  meter  and  300  meters  is  significant, 
but  the  difference  between  1  meter  and  5  meters  is  not. 

An  equation  for  terminal  velocity  could  be  built  into  the  data  entry  program  that 
takes  the  weight  of  the  chosen  mine  into  account  and  asks  the  user  if  the  water  depth  is 
less  than  that  which  would  produce  terminal  velocity  in  the  water  phase.  Assuming  no 
rotation  rate  was  chosen,  if  the  depth  was  known  to  be  greater  than  required  for  terminal 
velocity,  a  depth  need  not  be  entered.  This  simplification  would  also  allow  the  user  to 
have  some  confidence  in  the  result  as  he  moves  about  the  area,  regardless  of  water  depth 
changes. 

The  model  could  be  revised  to  provide  a  range  of  values  for  burial  depth,  based  on 
a  initial  attitude  of  90°  and  0°.  In  this  way,  the  uncertainty  of  initial  attitude  and  rotation 
rate  would  be  eliminated  and  a  more  realistic  range  of  values  would  be  produced  This 
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may  seem  like  a  reduction  in  sophistication  of  the  numerical  model,  but  the  reality  is  that 
the  exact  burial  depth  will  never  be  known  due  to  the  unpredictability  of  the  attitude  of 
the  mine  as  it  encounters  the  sediment  interface.  This  would  also  eliminate  any  effects 
due  to  currents  or  winds,  since  the  primary  effect  of  these  influences  would  be  on  the 
attitude  of  the  mine. 

Further  investigation  is  warranted  to  simplify  the  sediment  profile  data  entry. 
Assuming  the  values  for  density  and  shear  strength  are  either  known  or  can  be  measured 
at  the  interface,  a  set  of  equations  should  be  derived  and  refined  to  create  the  remainder 
of  the  profile.  If  this  were  an  option  in  the  model,  while  still  allowing  the  user  to  enter 
the  entire  profile  if  known,  it  would  substantially  increase  the  usefulness  and  precision  of 
the  model.  The  few  cases  discussed  here  and  the  equations  derived  from  that  limited  data 
set  are  encouraging,  and  may  indicate  that  such  equations  are  possible  and  beneficial. 

A  more  precise  method  for  simplifying  the  sediment  data  ingest  would  be  to 
obtain  samples  from  areas  where  mine  clearing  operations  are  likely  to  be  accomplished. 
The  relationship  between  density  and  shear  strength  for  each  sample  would  then  be 
computed  using  the  a  and  p  coefficients  discussed  previously.  The  user  could  determine 
sediment  type  and  density  using  available  environmental  atlases  and  the  model  would 
synthesize  the  profile. 
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